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Cellular systems are generally robust against fluctuations of intracellular parameters such as gene 
expression level. However, little is known about expression limits of genes required to halt cellular 
systems. In this study, using the fission yeast Schizosaccharomyces pombe, we developed a genetic 
'tug-of-war' (gTOW) method to assess the overexpression limit of certain genes. Using gTOW, we 
determined copy number limits for 31 cell-cycle regulators; the limits varied from 1 to >100. 
Comparison with orthologs of the budding yeast Saccharomyces cerevisiae suggested the presence of 
a conserved fragile core in the eukaryotic cell cycle. Robustness profiles of networks regulating 
cytokinesis in both yeasts (septation-initiation network (SIN) and mitotic exit network (MEN)) were 
quite different, probably reflecting differences in their physiologic functions. Fragility in the 
regulation of GTPase spgl was due to dosage imbalance against GTPase-activating protein (GAP) 
byr4. Using the gTOW data, we modified a mathematical model and successfully reproduced the 
robustness of the S. pombe cell cycle with the model. 
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Introduction 

Intracellular parameters such as gene expression require 
optimization, such that cellular functions may be performed 
effectively (Alon et al 1999; Zaslaver et al 2004; Dekel and 
Alon, 2005; Wagner, 2005). Fluctuations in these parameters 
lead to various cellular defects. Overexpression of genes 
involved in proliferation of cancer cells due to gene amplifica- 
tion is a prime example (Albertson, 2006) . On the other hand, 
in order to maintain cellular functions despite environmental 
change, mutation, and noise in intracellular biochemical 
reactions, these parameters may have certain permissible 
ranges, a characteristic termed robustness, which is com- 
monly observed in various cellular systems (Barkai and 
Leibler, 1997; Little et al 1999; von Dassow et al 2000). 

We previously reported a method designated genetic 'tug-of- 
war (gTOW),' by which we can measure the limit of gene 
overexpression in the budding yeast Saccharomyces cerevisiae 
(Moriya et al 2006). In gTOW, a target gene with its native 
regulatory regions is cloned into a particular plasmid, and the 
copy number of the plasmid is increased by genetic selection. 
Next, the copy number is measured just before the cellular 
system halts (i.e., the cell dies), such that the overexpression 
limit of the target gene is evaluated as the gene copy number 
limit. As the gene copy number increases, relative over- 
expression of the gene is expected. If we can measure the copy 
number limit of gene overexpression, we can evaluate the 
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degree to which the cellular system resists overexpression of 
the target gene (namely, robustness versus gene overexpres- 
sion) . 

Using gTOW, we previously measured the copy number 
limits of 30 cell-cycle regulators in budding yeast (Moriya et al 
2006) . The data were used to reveal the robustness profile of 
the cell-cycle regulatory system, and to evaluate and refine the 
integrative mathematical model of the budding yeast cell cycle 
(Moriya et al 2006; Kaizu et al 2010). The fission yeast 
Schizosaccharomyces pombe is distantly related to S. cerevisiae 
(Sipiczki, 2000), and like S. cerevisiae, is an established model 
eukaryote for the study of the molecular biology of the cell 
cycle (Egel, 2004) . In this study, we developed a gTOW method 
using S. pombe and determined the copy number limits of 31 
cell-cycle regulators. The data thus obtained were used to 
compare the robustness profiles of budding and fission yeasts in 
order to reveal the conserved/non-conserved properties of the 
eukaryotic cell cycle. The data were also used to constitute an 
integrative mathematical model of the fission yeast cell cycle. 

Results 

Development of plasmid vectors used for 
S. pombe gTOW 

The scheme of S. pombe gTOW is described in Figure 1A. The 
plasmid for gTOW must have the following properties: (1) the 
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copy number in each cell is multiple; (2) the copy number is 
diverse in each cell; and (3) the plasmid contains a gene with a 
selection bias for increasing the plasmid copy number [leu2d 
in the case of S. cerevisiae gTOW) . Because S. pombe possesses 
no plasmid with the above properties, we first constructed 
plasmid vectors for use in S. pombe gTOW. 

Given below is a brief summary of the construction of the 
plasmid (details are given in Supplementary information) . We 
used pDblet (Brun et al, 1995) for the backbone of the plasmid 
because it satisfies properties 1 and 2 above, and built S. 
cerevisiae LEU2 with promoter deletions into the plasmid to 
satisfy property 3. We first used S. cerevisiae leu2d (a LEU2 
allele used in S. cerevisiae gTOW). In S. cerevisiae, the selection 
bias of leu2d to increase the plasmid copy number can be 
controlled by the leucine concentration in the medium (Moriya 
et al, 2006). However, this proved impossible in S. pombe 
because the final copy numbers of the leu2d plasmid 
(pTOWsp-H) were the same in media containing different 
leucine concentrations (data not shown) . This is probably due 
to differences in the leucine-sensing mechanism between the 



yeasts. We thus optimized lengths of the LEU2 promoters 
(Supplementary Table S4) and developed three independent 
gTOW vectors (pTOWsp-L, pTOWsp-M, and pTOWsp-H) with 
different biases to increase plasmid copy numbers in leulA 
S. pombe cells (Supplementary Figure SI). Average plasmid 
copy numbers of the plasmids in the leal A cell in the leucine + 
media (uracil- selection) were 8.5 ±3.3 (all three plasmids 
showed no difference), and in the leucine— media were 
-41.3 ±13.0, 55.4 ±11.3, and 116.8 ±7.0, respectively 
(Figure IB; Supplementary Table S4). Differences between 
the gTOW vectors in S. cerevisiae and S. pombe are shown in 
Supplementary Table S5. 

In each plasmid, Ura4-EGFP was expressed, such that 
distribution of the plasmid copy number within populations 
could be monitored. Distribution of the cells for each of the 
plasmids analyzed by flow cytometry is shown in Figure 1C. 
In the leucine + condition (orange in Figure 1C), the plasmids 
tended to become lost because a portion of the cells showed 
the same fluorescence as the negative control (gray in 
Figure 1C). The average GFP fluorescence and copy number 
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Figure 1 gTOW in Schizosaccharomyces pombe. (A) The four steps in the gTOW experiment. (1 ) Each target gene is cloned into a plasmid for gTOW. (2) The plasmid 
is introduced into the leu1Aura4A cell (first selection is uracil). The plasmid copy number becomes ~ 10 per cell due to ars3002 x 2 ORI. If the copy number limit of the 
target gene is < 10 per cell, then the plasmid copy number is determined according to the limit. (3) The cells are transferred into leucine- media. In this condition, cells 
with a high plasmid copy number are selected according to the strength of expression of the leucine synthesizing enzyme gene LEU2. On the other hand, if the increase 
in the target gene (i.e., overexpression) is toxic to the cell, then the plasmid copy number must be less than the toxicity limit. The resulting tug-of-war of both selection 
biases determines the plasmid copy number within the cell, the number being close to the copy number limit of overexpression of the target gene. (4) To estimate the 
copy number limit of the target gene, the growth rate, GFP fluorescence of the cell, and the plasmid copy number within the cell are measured. In gTOW, the growth rate 
of the cell with the plasmid and the plasmid copy number within the cell should be mutually related (Moriya et al, 2006). We observed this trend in S. pombe, but 
measuring growth rate with a microplate reader, as described previously (Moriya etal, 2006), was unstable in this species. We thus cannot show the growth rate of the 
cell. (B) Comparison of the plasmid copy number determined by real-time PCR with mean GFP fluorescence of the cell determined by flow cytometry. A linear correlation 
between the plasmid copy number and the GFP fluorescence is observed. (C) Flow cytometric analysis of cellular distribution for each pTOW plasmid. pHM002 
represents the control without GFP. In the leucine + condition, no difference is observed among pTOWsp-L, -M, and -H. 
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showed good correlation (Figure IB), which indicates that GFP 
fluorescence is also useful in estimating the plasmid copy 
number within the cell. 



Measuring copy number limits of cell division 
cycle regulator (cdc) genes using gTOW 

We first selected 31 cdc genes (except byr4 + , Table I) and 
measured their copy number limits using gTOW. These genes 
contain orthologs of the budding yeast cdc genes, which we 
had previously analyzed (Moriya et al, 2006). They also 
contained genes implemented into an integrative mathema- 
tical model described below. In gTOW, one needs to measure 
by how much the expression of a target gene can be increased 
from the native expression level; hence, whole regulatory 
elements on the 5' and 3' regions of the target ORF should be 
cloned into the plasmids. We thus cloned each target gene just 
before and after the neighboring ORFs (the case of cdc2 is 
shown in Supplementary Figure S3 as an example). In 
addition, we introduced frameshift mutations in five genes 
just after their start codons to test whether proteins expressed 
by these genes determine copy number limits (the structure of 
the frameshift mutant of cdc2 is shown in Supplementary 
Figure S4 as an example) . 

Each gene was cloned into pTOWsp-L, -M, and -H vectors, 
and at least two independent clones were studied for each 
gene. S. pombe cells with each plasmid were grown in 



Table I Schizosacchawmyes pombe genes analyzed in this study 



No. 


Gene 


Description 


1 


arkl + 


Aurora kinase 


2 


cdc2 + 


Cyclin-dependent kinase 


3 


cdc7 + 


Septation kinase 


4 


cdclO + 


Gl phase transcription factor 


5 


cdcl3 + 


B-type cyclin, M phase 


6 


cdcl6 + 


With Byr4, a two-component GEF for the GTPase Spgl 


7 


cdcl8 + 


S-phase initiator/MCM loader 


8 


cdc25 + 


CDK tyrosine phosphatase 


9 


chkl + 


Checkpoint kinase, damage response 


10 


cigl + 


B-type cyclin, function not clear 


11 


cig2 + 


B-type cyclin in S phase 


12 


clpl + 


Phosphatase, involved in septation 


13 


cskl + 


CDK-activating kinase 


14 


cutl + 


Promotes anaphase, Separase 


15 


cut2 + 


APC substrate, blocks Cutl, Securin 


16 


dfpl + 


Regulatory subunit of CDC7-type kinase 


17 


fkh2 + 


Fork head transcription factor 


18 


hskl + 


DBF4-dependent Kinase 


19 


mikl + 


CDK tyrosine kinase, overlaps with weel + 


20 


plol + 


Polo-like kinase 


21 


pucl + 


Cyclin, acts in Gl at cell-cycle entry/ exit 


22 


pyp3 + 


CDK tyrosine phosphatase 


23 


rasl + 


Ras homolog 


24 


resl + 


Partner to CdclO transcription factor 


25 


res2 + 


Partner to CdclO transcription factor 


26 


rural + 


CDK inhibitor, functionally similar to Sicl 


27 


sid2 + 


Ser/Thr kinase, functions as part of a network in 






mitosis 


28 


slpl + 


APC activator 


29 


spgl + 


GTPase, involved in septation and mitotic exit 


30 


srwl + 


APC activator, WD protein 


31 


weel + 


CDK tyrosine kinase 


32 


byr4 + 


With Cdcl6p, a two-component GEF for the 






GTPase Spgl 



leucine + and leucine— (uracil—) media. After 48 h of 
cultivation for pTOWsp-L and -M plasmids and 72 h for 
pTOWsp-H plasmid, plasmid copy numbers within the cells 
and GFP fluorescence of the cells were measured and cellular 
morphologies were observed. Cell size distribution was also 
monitored as forward scattering (FSC) by flow cytometry, and 
some examples are shown in Figure 2. Using combinations of 
leucine + and leucine— conditions with three vectors, the 
copy number limit of each gene could be determined precisely 
(second and third columns in Figure 2). Using fluorescent 
microscopy, we could assess how cells with high target gene 
copy numbers (brighter GFP fluorescence) died (fifth and sixth 
columns in Figure 2). For example, in the case of cdc25, the 
copy number limit could be determined by the pTOWsp-M 
vector, and cells became smaller in the leucine— condition. In 
the case of rural, both cells and nuclei enlarged, while in the 
case of spgl, cells showed multiple septation. The pTOWsp-H 
vector for these genes showed a very strong bias for cells with 
the plasmid in the leucine— condition (red graphs in the FACS 
histogram in the third column of Figure 2) . 

In addition, when we analyzed the gTOW experiments by 
time-lapse microscopy, we could assess how cells with high 
target gene copy numbers died. Figure 3 shows an example of 
time4apse observation of a gTOW experiment with the target 
gene mikl + , which encodes a mitotic inhibitory kinase. 
Daughter cells with unequally distributed plasmid copy 
numbers had different fates; cells with low plasmid copy 
numbers barely grew, and after two cellular divisions growth 
ceased completely (black arrowhead in Figure 3). Cells with 
high plasmid copy numbers grew faster because they produced 
leucine but were unable to divide, resulting in elongation and 
eventual death (white arrowhead in Figure 3). On the other 
hand, cells with a balanced plasmid copy number maintained 
normal growth and division and accumulated in the media 
(yellow arrowhead in Figure 3). 

The copy number limits varied from 1 (one extra copy of the 
endogenous copy number kills the cell) to >100 (Figure 4), 
indicating that the cell-cycle regulatory system has a different 
robustness against overexpression of each gene, as observed in 
S. cerevisiae cell-cycle regulation (Moriya et al, 2006). The 
frameshift mutants generally showed much higher copy 
number limits than the wild types (orange bars in Figure 4), 
but the frameshift mutants of weel and ruml had lower limits 
than the pTOWsp-H vector alone, probably because their 
frameshifts do not completely shut down protein expression. 
The morphology of cells with high copy frameshift mutant 
genes was similar to that of cells with high copies of their wild 
types (Supplementary Figure S6) . 



Comparison of robustness profiles of cell-cycle 
regulation in S. cerevisiae and S. pombe 

We next compared the copy number limits of the cdc genes 
determined in this study with that of the S. cerevisiae homologs 
determined previously (Moriya et al, 2006; Figure 5). S. pombe 
genes generally had lower copy number limits (average=52) 
than those of S. cerevisiae (average=87). We previously 
reported that cdc genes with lower limits constituted a fragile 
core that directly regulates the activity of B-type cyclin Cdk 
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Figure 2 Examples of gTOW experiments with cdc genes in S. pombe. The experimental results obtained by gTOW with empty vector, cdc2+ , cdc25+ , ruml + , 
and spgl + are shown (enlarged part of figure is shown in Supplementary Figure S5). First column: The plasmid copy number (blue graph, right vertical axis) and the 
mean GFP fluorescence (green graph, left vertical axis) in gTOW in leucine + (indicated as ' + ') and leucine- with each vector (pTOWsp-L, -M, and -H) are shown 
(indicated as 'L-' 'M-,' and 'H-,' respectively). The original data are shown in Supplementary Tables S6 and S7. Second column: Cell distribution with GFP 
fluorescence in gTOW experiments in leucine + (orange graph) and leucine- conditions with -L vector (green graph), -M vector (blue graph), and -H vector (red 
graph) are shown. Some experiments showed bimodality distributions. This is probably due to the plasmid loss, which is a property of the ars-based plasmid. Third 
column: Scatter plot between GFP fluorescence (FL-1) and cell size (FSC). Cells were cultured in the leucine + condition. Fourth column: Fluorescent microscopic 
image of the cells cultured in the leucine + condition. Nucleus and septum are stained. Fifth column: GFP fluorescence of the cells in the fourth column. The brightness 
reflects the plasmid copy number within the cell. The whole data set of above analysis for the genes analyzed in this study is provided upon request. 



(Cyclin-dependent kinase) complexes (Moriya et al, 2006). 
The orthologs that constitute the core structure of S. pombe 
(Cdcl3, Ruml, and Weel) also demonstrated low limits, 
which may indicate the presence of a conserved fragile core in 
eukaryotic cell-cycle regulation. In contrast, the orthologous 
components of the mitotic exit network (MEN) of S. cerevisiae 
showed very different limits from those in the septation- 
initiation network (SIN) of S. pombe, although both networks 
have conserved architecture (Bardin and Amon, 2001; Krapp 
and Simanis, 2008) . This might reflect the different physiologic 

4 Molecular Systems Biology 2011 



functions in the two yeasts (namely, budding and fission), 
although these components have the same origin. These 
results suggest that comparison of robustness profiles is a 
useful way to reveal the conserved/non-conserved properties 
of cellular systems. 

Dosage imbalance causes very low limit of spgl 

We recently demonstrated that the M-phase phosphatase gene 
CDC14 has a very low limit is because of the dosage imbalance 
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Figure 3 Time-lapse observation of gTOW experiment with the target gene 
mik1 + . (A) FY7652 cells harboring the pTOWsp-H plasmid with the mitotic CDK 
inhibitor mik1 + were cultivated in EMM without leucine and were observed 
under fluorescent microscopy. The plasmid copy number in each cell was 
estimated from the intensity of GFP fluorescence. The black arrowhead with a 
white frame indicates a cell with very low plasmid copy number to synthesize a 
sufficient amount of leucine to support growth. The cell stops proliferating after 
the second cell division (frame 6). The white arrowhead indicates a cell with a 
sufficiently high plasmid copy number to synthesize enough leucine to support 
rapid cellular growth, but this number is beyond the copy number limit of mik1 + . 
The cell thus cannot divide, becomes elongated (frames 3-5), and eventually 
dies (frame 6). The yellow arrowhead indicates cells with balanced plasmid copy 
number; the cells can grow and divide because they synthesize sufficient leucine. 
Cells with a plasmid copy number close to the limit of mik1 + thus become 
concentrated during cultivation. A video clip is provided in Supplementary 
information. (B) Assumed cellular situations in gTOW experiment are shown (A). 
The arrowheads correspond to the arrowheads in (A). 
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Figure 4 Copy number limits of 30 cell-cycle regulators in S. pombe. The 
green bar indicates the plasmid copy number obtained in the gTOW experiment 
with each vector. The blue bar indicates the maximum plasmid copy number 
for each target gene obtained in the gTOW experiment with three different 
vectors. The orange bar indicates the maximum plasmid copy number for each 
target gene with frameshift mutation. The original data are shown in 
Supplementary Table S6. Because the cell has an endogenous copy 
of the target gene, the copy number limit for the target gene is the plasmid 
copy number plus 1. 



between CDC14 and its inhibitor gene, NET1 (Kaizu et al, 
2010). Hence, we investigated whether the low limit of spgl, 
the lowest copy number gene in S. pombe cell-cycle regulators 
evaluated in this study, was also due to dosage imbalance, spgl 
encodes a small GTPase involved in the initiation of cellular 
septation; this activity is regulated by the bipartite GTPase- 
activating protein (GAP), which is encoded by both cdcl6 and 
byr4 (Furge et al, 1998; Figure 6A). The copy number limit of 
cdcl6 is quite high (Figure 4), and cells with high copy cdcl6 
did not show any obvious phenotype (data not shown). In 
contrast, the limit of byr4 was very low (1 .8 ± 0.2, measured by 
pTOWsp-L in the leucine + condition) . And cellular lethality 
due to overexpression of spgl is nullified by the simultaneous 
overexpression of byr4 (Furge et al, 1998). We thus assume 
that proper activation of Spgl is in delicate balance with Byr4, 
and hence, the limits of both spgl and byr4 are quite low. 
To investigate whether spgl and byr4 were in dosage balance, 
we performed two-dimensional (2D) gTOW (Kaizu et al, 2010). 
As shown in Figure 6B, byr4 limits increased only when 
multicopy spgl was supplied by another plasmid. Moreover, 
their copy numbers were balanced (Figure 6B, dotted line). 
Microscopic observation also confirmed our assumption that 
cells can divide normally only when spgl and byr4 expressions 
are balanced (Figure 6C) . These results strongly support the 
theory that spgl has a very low limit because of the sensitive 
balance required for GTPase activity. 



Development of an integrative mathematical model 
that reproduces the robustness of cell-cycle 
regulation in S. pombe 

Robustness can be a measure of plausibility in mathematical 
models of biological networks (Morohashi et al, 2002). We 
previously used copy number limits obtained for S. cerevisiae 
gTOW to evaluate and refine an integrative mathematical 
model of the S. cerevisiae cell cycle (Chen et al, 2004; Moriya 
et al, 2006; Kaizu et al, 2010). We thus evaluated and refined 
the integrative mathematical model, such that the model 
reproduces robustness (copy number limits of genes) of the S. 
pombe cell cycle obtained from this study. We previously 
published a mathematical model for S. pombe cell-cycle 
regulation (Novak and Tyson, 1997; Sveiczer et al, 2000), 
and have independently modified the model from the gTOW 
experiment, such that the model reproduces more published 
experimental results. Here, we designated the model 'basic 
model,' the whole structure of which is shown in Figure 7 A 
(green components; the simplified structure is shown in 
Supplementary Figure S8) . We first investigated to what degree 
the core model reproduced the copy number limits obtained in 
this study (see Supplementary information for details). 
Although we did not use the gTOW data for development of 
the model, the model predicted the experimental data well 
(Figure 7B, green circles), and appeared to capture the 
robustness of fission yeast cell-cycle regulation. To describe 
the gTOW data more extensively, we modified the basic model 
by adding some important regulators (M-phase phosphatase 
Clpl, cyclins Cigl and Pucl) and regulations, and the 
parameters were optimized by-hand parameter adjustments 
(we designated the model 'gTOW model,' shown in red in 
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Figure 5 Comparison of robustness profiles of cell-cycle regulation between S. cerevisiae and S. pombe. The dotted blue line connects the functional orthologs within 
the 'fragile core' conserved in both yeasts (left). The dotted orange line connects the functional orthologs involved in the cytokinesis regulatory pathway of both yeasts 
(MEN and SIN). Some components involved in both MEN and SIN are omitted from the diagrams because their copy number limits were not measured. Data for 
S. cerevisiae are obtained from a previous study (Moriya et al, 2006). The original data are shown in Supplementary Table S8. 



Figure 7 A; see also Supplementary Figure S8). The model 
successfully reproduced the copy number limits obtained in 
this study (Figure 7B, orange squares). In addition, the model 
could also reproduce 42 already published mutant behaviors 
(Supplementary Table S10). The result of time-course simula- 
tion of the model is shown in Figure 7C. As a result, we can 
claim that the presented 'gTOW' model is the most detailed 
model of fission yeast cell-cycle regulation so far. 



Increase of each Cdc2/cyclin complex, but not the 
depletion of Cdc2/Cdc13 through the competition 
for Cdc2, determines the upper limits of cyclins 
(cigl, cig2, and pud) 

S. pombe contains four cyclins (Cdcl3, Cigl, Cig2, and Pucl), 
all of which bind to Cdc2 to form a Cdc2/cyclin complex that 
performs regulatory functions on Cdc2 (Figure 8A). Cdk's (like 
Cdc2) are thought to be in excess of cyclins, but in S. cerevisiae 
it is not tremendously so (Cross et al, 2002). It is thus possible 
that if cyclins compete with each other for Cdc2, over- 
expression of one of the cyclins would interfere in formation 
of other Cdc2/cyclin complexes. Especially since deletion of 
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Cdcl3 can lead to lethality, there is a chance that over- 
expression of the other cyclins is deleterious for the cells 
because Cdc2 is titrated away from Cdcl3 (Figure 8A-2, case2) . 
In other words, the overexpression limit of each cyclin may be 
determined by the increased Cdc2/cyclin activity itself (Figure 
8A-2, easel) or by the reduction of Cdc2/Cdcl3 complex. 
To answer this question, we first used the mathematical model 
to investigate which conditions involve cyclins in competitive 
situations, and the experimental means to test these condi- 
tions. In the gTOW model (and in the basic model), each cyclin 
is considered to form a Cdc2/cyclin complex immediately after 
its expression; thus, no monomer cyclins or Cdc2 would be 
available and one cannot simulate competition with these 
models. Hence, we further modified the gTOW model, such 
that the behavior of monomer cyclins and Cdc2 could be 
simulated (green components in Supplementary Figure S8). 

We first investigated the parameters of the condition in 
which cyclin competition occurs using this model, and 
observed that the condition is determined by a combination 
of the dissociation constant of Cdc2/cyclin and total Cdc2 
level, and that this condition can be experimentally evaluated 
by measuring the copy number limit of each cyclin in the Cdc2 
overexpression condition (Supplementary Table S12; Figure 
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8A-3). If the limit is decreased, then the limit should be 
determined by an increase in the Cdc2/cyclin activity itself; 
and if the limit is increased, then the limit should be 
determined by the reduction of Cdc2/Cdcl3 (competition). 
Figure 8B shows the copy number limits of cigl , dg2, and pucl 
with multicopy cdc2; in all cases, the limits were decreased. We 
thus concluded that upper limits of cyclin genes are 
determined because of the high Cdc2/cyclin (Cigl, Cig2, and 
Pucl) activity itself, but not the depletion of Cdc2/Cdcl3. This 
conclusion is supported by the phenotypic observation of 
gTOW-agi + , -cig2 + , and pucl + cells (Supplementary 
Figure S7). The cells did not show any elongated cellular 
morphology, which is a typical cellular phenotype when 
Cdcl3-Cdc2 activity is depleted. 



Comparison of 2D-gTOW data in silico and in vivo 
revealed uncovered regulations in the cell cycle 

One of the practical values of the existence of the integrative 
mathematical model, just as we developed in this study, is that 
we can define the limit of our understanding and find novel 
biological knowledge using the model predictions as refer- 
ences of experimental results. We thus measured the copy 
number limits of the cell-cycle regulators in the deletion 
strains, and compared the data with the model prediction in 
Supplementary Table Sll. This £ 2D-gTOW experiment' was 
performed with the deletion strain of each of four genes {cigl, 
dpi, pucl, and srwl). We chose these genes because the 
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deletion of them does not give so much impact for the upper 
limits of the other cell-cycle regulators in the model (Figure 9A; 
Supplementary Table Sll). We thus expected that we could 
reveal uncovered regulations. 

In fact, we found some discrepancies between the gTOW 
experimental data and the gTOW model predictions 
(Figure 9A). The rural limit in cigl A, the slpl limit in dpi A, 
and cdc25 limit in srwl A were much lower than the model 
prediction. The pyp3 limit in srwl A was much higher than the 
model prediction. These results suggested additional regula- 
tions and wrong assumptions within the current gTOW model 
(summarized in Figure 9B) . 

Discussion 
gTOW in S. pombe 

In this study, we developed a gTOW method in the fission 
yeast S. pombe. gTOW vectors in S. pombe have different 
biases for increasing the plasmid copy number. These vectors 
also have GFP, such that distribution of the plasmid copy 
numbers can be monitored by flow cytometry, or the plasmid 
copy number in a single cell can be estimated under 
fluorescent microscopy. In S. pombe, the nmtl promoter has 
been often used to overexpress genes, but the upper limit is 
difficult to determine by this method. Thus, gTOW provides 
another means for overexpression of genes and for obtaining 
unique data to investigate the robustness of cellular systems in 
S. pombe. 

Molecular Systems Biology 2011 7 



Limits of fission yeast cell-cycle genes 
H Moriya et al 




B 



cdc10+ 



wee1+ y 



srwl 



cdc18+ 



slp1+ 



rumU 



pyp3+ 




cdc25+ 



mik1+ 



V gTOW experiment 
# DNA replication model 
■— gTOW model 




200 



Figure 7 Mathematical model reproducing gTOW data. (A) Whole structure of the mathematical model of the fission yeast cell cycle developed in this study (gTOW 
model) given with Systems Biology Graphical Notation (SBGN) using CellDesigner 4.1 . Pink colored components are the ones added to the 'basic model' to make the 
'gTOW model'. CellDesigner file and Systems Biology Markup Language (SBML) file are provided in Supplementary information. (B) Comparison of the copy number 
limits of cell-cycle regulators between the data obtained by gTOW and prediction of the mathematical model. Scale of the axis=log2. (C) Time-course simulation result of 
the gTOW model. 



Robustness and fragility in cell-cycle regulation in 
S. pombe and S. cerevisiae 

S. cerevisiae and S. pombe appear to have conserved fragility in 
the subsystem that directly regulates B-type cyclin Cdkl 
activity (Figure 5) . Because these species are distantly related, 
with an evolutional distance of several hundred million years 
(Sipiczki, 2000), fragility in B-type Cdk regulation would be 
conserved among eukaryotes. In contrast, networks regulating 
cytokinesis (MEN in S. cerevisiae and SIN in S. pombe) showed 
different robustness profiles (Figure 5), which might reflect 
their different systems of cellular division (budding and 



fission). Hence, comparison of robustness profiles would be 
useful to reveal the conserved properties of cellular systems. 

In the cell-cycle regulation of S. cerevisiae, the Cdkl- 
counteracting phosphatase gene CDC14 had the lowest copy 
number limit, the reason being dosage imbalance against the 
inhibitor gene NET1 (Kaizu et al, 2010). The relatively high 
copy number limit of dpi, the S. pombe ortholog of CDC14, 
may be explained by the fact that Clpl does not have a 
stoichiometric inhibitor (it seems there is an ortholog of Netl, 
but acts differently; Jin et al, 2007) . In this study, we showed 
that the small GTPase gene spgl , the gene with the lowest copy 
number limit among the S. pombe cell-cycle regulators, was 
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also in dosage balance with the GAP gene byr4 (Figure 6). 
We therefore believe that dosage imbalance is a general 
mechanism for inducing cellular fragility in case of gene 
overexpression. 

Integrative mathematical model of cell-cycle 
regulation in S. pombe 

In this study, we evaluated and modified mathematical models 
constituted from molecular knowledge using gTOW data. 
Although the initial model (the basic model) was indepen- 
dently developed from gTOW experiments, the model pre- 
dicted the experimental results well (Figure 7 A) . This model is 
a quantitative model in which biochemical parameters had not 
been experimentally determined, but it already appears to 
have captured the robustness characteristics of the S. pombe 
cell cycle. 

By addition of certain components to the model and 
optimization of parameters, we successfully developed a 
model (the gTOW model) to reproduce the gTOW data 
obtained in this study (Figure 7 A) and the behavior of 42 
other mutants (Supplementary Table S10). The model 
provided us with two types of predictions that should be 
evaluated by further research. One of these predictions 
concerns the effect of gene overexpression on phenotype. 
With simulations using the gTOW model, three different cell- 
cycle defects were observed due to overexpression of each 
gene (Supplementary Figure Sll). Interestingly, the model 
predicted that ruml overexpression results in three indepen- 
dent phenotypes according to the level of overexpression 
(Supplementary Figure SI 2). 
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The other prediction concerns the copy number limits of 
each gene in its deletion strains (Supplementary Table Sll). 
Upon comparison with the data in vivo, this type of prediction 
is quite useful to find additional regulations that are not 
implemented in the model (Kaizu et al, 2010) . In this study, we 
measured upper limits of cell-cycle regulators in four deletion 
mutant strains, which suggest some novel regulations in the 
network lacking in our model (Figure 9) . We thus showed that 
the prediction of the model was quite useful reference to 
indicate the existence of additional (unknown) regulations 
conferring the robustness of the cell cycle. We propose that the 
following potential regulatory mechanisms provide an 
explanation for the 2D-gTOW data in Figure 9. 

(1) Low ruml limit in cigl A: Cigl might be more important for 
Ruml degradation than it is assumed in the model. 

(2) Low slpl limit in clplA: Clpl might negatively regulate 
Slpl activity or the level that is absent in clpl A back- 
ground. 

(3) Low cdc25 in srwlA: Srwl might be important for 
degradation of Cdc25. 

(4) Low pyp3 limit in srwl A in the model but not in the 
experiment: Pyp3 might be less important as a Cdc25 
backup phosphatase than it is assumed in the model. 

At this point, we did not modify the gTOW model with these 
suggested regulatory interactions because they are still 
hypothetical. These hypothetical regulations should be tested 
with additional experiments and the model will be refined with 
the regulations, just as we did in the budding yeast case (Kaizu 
et al 2010). 
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Evaluation and refinement of integrative 
mathematical model using robustness information 

A considerable body of knowledge has been accumulated for 
the fission yeast S. pombe and for the molecular details of the 
regulation of its cell cycle (Nurse, 2002), and this knowledge 
has been integrated into mathematical models before (Novak 
and Tyson, 1997; Sveiczer et al, 2000), and further developed 
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in this study. These integrative mathematical models are 
expected to have very important roles in future life sciences 
because they can be used to predict behaviors of complicated 
biological systems that are difficult to handle intuitively. 

These models are generally constructed according to the 
accumulated knowledge on the interactions between relevant 
molecules, but information on biochemical parameters needed 
to describe the strength and speed of interactions is usually 
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lacking. Although these parameters are important to guarantee 
the plausibility of these models, it is sometimes technically 
difficult to obtain these biochemical parameters experimentally. 

On the other hand, it is considered that these parameters 
should have some permissible ranges because biological 
systems have robustness (Barkai and Leibler, 1997; Little 
et al, 1999; von Dassow et al, 2000). Because the permissible 
ranges have emerged from interactions of components within 
the system, they are useful indicators for evaluation of the 
network structures of biological models (Morohashi et al, 
2002; Moriya et al, 2006; Kaizu et al, 2010). In this study, we 
performed experiments to measure the permissible range of 
gene overexpression and used the data obtained to evaluate 
and modify the mathematical model. We believe this is a useful 
scheme for development of integrative mathematical models 
in which direct experimental measurement of biochemical 
parameters is difficult. 

Materials and methods 

Strains 

Saccharomyces cerevisiae strain BY4741 (MATa, his3Al, leu2A0, 
metl5A0, ura3A0) (Brachmann et al, 1998) was used for plasmid 
construction by the gap-repair cloning method (Oldenburg et al, 1997) . 
Schizosaccharomyces pombe strain FY7652 [h- leul-32 ura4-D18) 
obtained from NBRP yeast was used for gTOW, apart from 2D-gTOW 
{byr4 versus spgl and cyclins versus cdc2). In 2D-gTOWof S. pombe, 
Sp286h+ (7i+ ade6 ura4-D18 leul-32), a haploid progeny of Sp286 
(BG_0000, Bioneer) was used. For clplA, ciglA, puclA, and srwlA 
strains, haploid progenies of BG_0465, BG_4780, BG_4172, and 
BG_0289 (Bioneer) were used, respectively. The Escherichia coli strain 
XLl-blue (Stratagene) was used for plasmid amplification. 



Media and growth conditions 

YPD and SC medium were used for cultivation of S. cerevisiae and LB 
broth for E. coli. S. pombe cells were cultured in a manner similar to 
that described by Moreno et al (1991). EMM medium and plates were 
prepared using EMM broth powder (MP Biomedicals) . Yeast transfor- 
mation was performed using the LiOAc method (Amberg et al, 2005). 
Yeasts and E. coli were cultivated at 30 and 37°C, respectively. 



Plasmids used in this study 

Plasmid vectors used in this study are listed in Supplementary Table 
SI. Construction details of the plasmids used in this study are 
described in Supplementary Information. The vectors pTOWsp-L, 
pTOWsp-M, pTOWsp-H, and pA6R are available from the National 
BioResource Project (http://yeast.lab.nig.ac.jp/nig/index_en.html) . 
The DNA sequences are provided in Supplementary information. 



Measurement of plasmid copy number in S. pombe 
cells 

S. pombe cells with pTOWsp containing target genes were grown in 
EMM without uracil, transferred into EMM without leucine and uracil, 
and then cultivated for 48 h (for cells with pTOWsp-L- and -M-based 
plasmids) or 72 h (for cells with pTOWsp-H-based plasmids). In 2D- 
gTOW shown in Figures 6B and 8B, pTOW and pA6R derivatives were 
introduced into Sp286h+ cells, and the cells were grown in EMM 
without uracil and adenine; they were then transferred into EMM 
without leucine, uracil, and adenine and cultivated for 48 h. EMM was 
prepared using EMM broth powder (MP Biomedicals). The plasmid 
dosage in a yeast cell was determined using real-time PCR, as 
described previously (Moriya et al, 2006), except that the pombe-1 



primer set (5'-CCCTCAGTTGCTTCTCCTAA-3' and 5'-GTGACTTCGTGA 
ATCAAGTG-3'), the ura4-3 primer set (5' -TCTTTGGCTACTGGTTCCTA-3' 
and 5'-TATGTAGTCGCTTTGAAGGT-3'), and the ade6-l primer set 
(5'-GCCCATCGCTTAAACATCA-3' and 5'-TGCATCGGGGTCAGTAAAT-3') 
were used for quantification of genomic DNA, pTOW plasmids, and 
pA6R plasmids, respectively. Data are averages of at least two 
independent experiments. 

The copy number determined by PCR using the ura4-3 primer is 
directly associated with the plasmid copy number in the cell because 
no sequence on the genomes of FY7652 or Sp286h + is amplified by 
the ura4-3 primer set. In contrast, the copy number determined by PCR 
using the ade6-l primer is one copy more than the plasmid copy 
number in the cell because an ade6 sequence on the genome of 
Sp286h + is amplified by the primer set. The plasmid copy number is 
shown in all figures and tables. The copy number of a target gene in the 
wild-type cell is one copy more than the plasmid copy number as a 
result of one endogenous copy on the genome. 



Analysis of S. pombe cells by flow cytometry 

Cells expressing EGFP were directly analyzed by flow cytometry 
(FACScalibur; Becton, Dickinson and Company). EGFP fluorescence 
was detected as FL-1, and cell size was monitored as FSC. The data 
were analyzed using the FlowJo7 software (http://www.flowjo.com). 



Microscopic observation of S. pombe cells 

Cells were fixed by 70 % ethanol for > 10 min, pelleted, and resuspended 
in the VECTASHIELD mounting medium with DAPI (1/2 diluted with 
water, Vector Laboratories) and Fluorescent Brighter 28 (final concentra- 
tion 10 |ig/ml; Sigma) . The cell suspension (2 ul) was spotted onto a glass 
slide and covered with an 18-mm coverslip. The coverslip was sealed 
with nail polish, and the slide was observed using the Leica DMI 6000B 
microscope (Leica Microsystems). GFP fluorescence, DAPI and Fluor- 
escent Brighter 28 fluorescence, and RFP fluorescence, were observed 
with GFP, A, and RFP filter cubes (Leica Microsystems), respectively. 



Computation 

Development of mathematical models and their simulations results are 
described in Supplementary information. The model files for XPP-AUT 
and CellDesigner 4.2 are attached in Supplementary information. The 
SBML file of the gTOW model shown in Figure 7A can be obtained from 
BioModels database (http://www.ebi.ac.uk/biomodels-main/) (ID: 
MODEL1111040000). 



Supplementary information 

Supplementary information is available at the Molecular Systems 
Biology website (www.nature.com/msb). 
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